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Abstract 

A complete classification is presented of quantum and classical super- 
integrable systems in E2 that allow tlie separation of variables in polar 
coordinates and admit an additional integral of motion of order three in 
the momentum. New quantum superintegrable systems are discovered for 
which the potential is expressed in terms of the sixth Painleve transcen- 
dent or in terms of the Weierstrass elliptic function. 



1 Introduction 

The purpose of this article is to obtain and classify all classical and quantum 
Hamiltonians H that allow the separation of variables in polar coordinates and 
admit a third order integral of motion Y . The system under study is character- 
ized by three conserved quantities: 

H^Pl±Pi + v{r,e) (1.1) 

X^Ll + 2S{e) (1.2) 

r= A,jk{Ll„p{pl} + {gi{x,y),pi} + {g2{x,y),p2} (1.3) 

where 

V{r,e) = R{r) + ^ (1.4) 

Here R{r) and S{9) are arbitrary functions and Aijk are real constants. The 
polar coordinates are defined as usual: x = rcosO and y = rsin9. 

In classical mechanics pi , p2 are the cartesian components of linear momen- 
tum and L3 is the two-dimensional angular momentum. In quantum mechanics 
we have 

d d d 

ox oy oO 

The curly brackets, {■,■}, in (II. 3p denote an anticommutator in quantum me- 
chanics. In classical mechanics we have {L\tPiP2} = '^L\PiP2- 
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This study is part of a general program devoted to superintegrable systems 
in classical and quantum mechanics. Roughly speaking, an integrable system 
is superintegrable if it allows more integrals of motion than degrees of freedom. 
For more precise definitions and an extensive bibliography see e.g. [D H]. A 
system in n dimensions is maximally superintegrable if it allows 2n — 1 inte- 
grals. In classical mechanics the integrals of motion must be well defined and 
functionally independent functions on phase space and typically at least one 
subset of n integrals (including the Hamiltonian) is in involution. In quantum 
mechanics the integrals should be well-defined linear operators in the enveloping 
algebra of the Heisenberg algebra with basis {xi,pi, h} for i = 1, ...,n and they 
should be algebraically independent (within the Jordan algebra generated by 
their anticommutators) . 

The majority of publications on superintegrability is devoted to the quadratic 
case when the integral of motion is quadratic in the momenta (see e.g. [3]-|10]). 
Quadratic superintegrability for one particle in a scalar potential is related to 
multiseparability in the Schrodinger equation or the Hamilton- Jacobi equation, 
in quantum or classical mechanics, respectively. 

More recently some of the interest has shifted to higher order integrability. 
An infinite family of superintegrable and exactly solvable systems in a Euclidean 
plane has been proposed [TT1[T2]. The potential depends on a real number k > 0. 
It has been conjectured that this sytem is superintegrable in quantum mechanics 
for all integer values of k with one integral of order 2 and the other of order 2k. 
For brevity we will call it the TTW model. The superintegrability of the TTW 
model has been so far confirmed for odd values of fc > 3 [13]. In the classical 
case all bounded trajectories are periodic [12] for all integer and rational values 
of k and superintegrability has been proven for such values of k [M]. Both the 
trajectories and the higher order integral of motion can be expressed in terms 
of Chebyshev polynomials. The generalisation of the TTW model in a three 
dimensional Euclidean space has been recently proposed |15j . 

A systematic study of integrable systems in classical mechanics with one 
third order integral of motion was initiated by Drach in a remarkable paper 
published in 1935 [TB]. He considered classical Hamiltonian mechanics in a 
two dimensional complex Euclidean plane and found 10 potentials allowing a 
third order integral. More recently it was shown that 7 of these systems are 
actually quadratically superintegrable and that the third order integral of mo- 
tion is a Poisson commutator of two independent second order ones [TT] [T8] . 
Quadratically integrable (and superintegrable) potentials coincide in classical 
and quantum mechanics. This is not necessarily the case when higher order 
integrals are involved [HI |20] . 

A systematic search for quantum and classical superintegrable systems in 
a real Euclidean plane with one third order integral of motion and one first 
order or second order one was started in [3TJ [22]. In [35] the second order 
integral of motion was chosen so as to assure separation of variables in cartesian 
coordinates. This lead to several new classical superintegrable systems but 
mainly to completely new quantum ones, in which the potential is expressed in 
terms of Painleve transcendents. The integrals of motion generate polynomial 
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algebras [5], [33]- [30]. Their representation theory was used to calculate the 
energy spectra and a relation with supersymmetric quantum mechanics was 
used to calculate wave functions [^-[30]. 

The existence of a third order integral of motion in quantum mechanics was 
investigated earlier [3T] and potentials expressed in terms of the Weierstrass 
function were obtained. 

In this article we continue with the classification of superintegrable systems 
and impose the conditions (1.1) to (1-4), i.e. separation of variables in polar 
coordinates. 

The determining equations for the existence of a third order integral are 
presented in Section 2 in polar coordinates. The possible form of the radial part 
of the potential is established in Section 3 and summed up in Theorem 1. The 
angular part is discussed in Section 4. Genuinely new superintegrable potentials 
are obtained when the radial part of the potential vanishes. The angular parts 
of the potential are expressed in terms of the sixth Painleve transcendent Pq or 
in terms of the Weierstrass elliptic function. The main results are summed up 
as theorems in Section 5. 



2 Determining equations of a third order inte- 
gral of motion in polar coordinates 

The quantities X and H commute in quantum mechanics and Poisson-commute 
in the classical case. The form of (|1.3p assures that all terms of order 4 and 3 
in the (Poisson) commutator [H, Y] = vanish. The vanishing of lower order 
terms provides the following determining equations for the functions gi , g2 and 
V in dlll]) and ([Ojl : 



If 

GiVr + G2Ve = — 



FiVrrr + F^VrrO + F^VrOO + FiVeee + rF^Vrr + {irFi 



-F.Wr 



2Fi 



-FaVee + (-F3 + 2Ci cos 9 + 2C2 sin d)Vr 



;F2 + 8D0 



{-2Ci sin e + 2C2 cost 



Ve 



(2.1) 



(Gi), 

{G2)0 
{G2)r 



^3FiVr + F2Ve 

= FsVr + SFiVe - ^ 



2{F2Vr + F3Ve) 



(Gi)e 



(2.2) 
(2.3) 

(2.4) 
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with: 



and 



Fi (6) = Ai cos 36* + A2 sin 36* + cos 61 + ^4 sin 0, (2.5) 

F2 (r, 0) = - (-3Ai sin 36* + 3A2 cos 36* - A3 sin 6 + A^ cos 9) 
r 

+ Bicos20 + B2sin26' + So, (2.6) 
i^3(r,0) = \{-iAi cos 36* -3^2 sin 36* + A3 cos 6* + sin 6*) 

+ -{-2Bi sin 261 + 2^2 cos 26*) 
r 

+ CiCos6i + C2sin6l, (2.7) 
F4 (r, 6*) 4r ( Ai sin 36* - A2 cos 36* - A3 sin 6* + yl4 cos 6*) 

+ ^(-Bi COS26'- B2 sin 26* + Bo) 

+ i(-Cisin6l + C2Cose') +£'0, (2.8) 



^ / a , • /3 n I a\ -51 sm 6^ + 32 cos / 
Gi (r, e*) = 51 cos + 52 sm 9, G2 (r, b*) = 



r 

The constants A^, B^, Q et Z?o are related to Aijk of (|1.3p : 

. ^030 — ^012 . Aq21 — A003 . 3ylo30 + A012 

^ 4 4 4 

4 3^003 + A021 A120 — A102 Am 

4 4 2 2 

Bq = ^^-^ —, Ci — v42io, (^2 = A201, -Do = A300 

The constants A^, B^, et Dq are real. Relations equivalent to (|2.1l) - (|2.4p were 
already obtained in [21] in cartesian coordinates. 

The third order constant of motion can be written in terms of A;, _Bi, Ci 
and Do as 

Y = DoLl + Ci{Llpi} + C2{LIp2} + BMpI+pD + B,{L3,pI~pI) 

+ 2B2{L3,pip2} + Aipi{pl - ipl) + A2P2{ipl - pI) + {A3P1 + AiP2){pl+pl) 

+ {9i[x,y),pi} + {g2{x,y),p2} (2.9) 

Under rotations the Hamiltonian (jl.ll) and the integral X (|1.2I) remain invari- 
ant but the third order integral Y (|1.3p or (|2.9p transforms into a new integral 
of the same form with new coefficients and new functions gi,g2- The constants 
Dq and Bq are singlets under rotations and hence invariants. The expressions 
{i?i,i?2}, {^1,^2} and {^3,^4} are doublets under rotations. 
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The existence of the third order integral Y depends on the compatibility of 
(|2.ip - (|2.4l) . Equation p.l|) establishes the difference between the quantum and 
the classical cases. We obtain the classical analog of (|2.1I) by setting ft n- 0. 

From (|2.2p - (|2.4[) we can deduce a compatibility condition for the potential 
V(r,9), namely the third order linear differential equation: 

= r^FaVrrr + (Sr^Fi - 2r^F2jVrr8 + (SFI - 2r^F3jVree + F2Vggg 

+ (2r^F3r + Gr^Fs - 2r^F2g - 3rFi)v;, + (2^20 - irF^ - 2r^F3r)Vee 
+ (6r^Fir + 187-3^4 - 2r^{F2r + i^ae) - SrFs + QFw)Vre 
+ (7-*i^3,r + er^Fs, + r2(6F3 - 2F2,e) - 4rF2e + iFwe)Vr 
+ (?,r^Firr + ISr^Fir + r2(18F4 - 2^3,9) - 7-(F2, + 4F3e) + F2gg^Ve (2.10) 

Compatibility of (|2.2l) - (|2.4p and (|2.ip imposes further conditions on V{r, 9), they 
are however nonlinear [21j. 

Our aim is to find all solutions of (|2.ip - (|2.4l) and thereby find all superin- 
tegrable classical and quantum systems that separate in polar coordinates and 
allow a third order integral Y. The existence of X (|1.2p is guaranteed by the 
form l|1.4p of the potential and is directly related to the separation of variables 
in polar coordinates. 

Specifying the radial dependance of (|2.6p to (|2.8p : 

F2ir,9) = -F2i{9)+F2oi0) (2.11) 
r 

F3{r,9) = ^F32{9) + ^Fsi{9) + F3oi9) (2.12) 

Fiir,9) = ^F^m + ^Fi2{e) + \Fii{9) + Z?o (2.13) 

we can deduce: 

Gi(r, 9) = 3Fi (r + -L5) - (|i + ^) 5 + m (2.14) 

where S{9) = -jp:S{9). The derivative with respect to r will be denoted by a 
d9 

prime so that R'{r) = —R{r). 

dr 

3 Radial term in the potential 

With the separation of the potential in polar coordinates (|1.4I) . (|2.10p can be 
expressed as: 

rri3i?(3) + m2R" + miR' + 7135^^) + 712^ + niS + n^S = (3.1) 
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for rrii = rrii (r, 9) and Ui = Ui (r, 9) in terms of the functions Fi (r, 0) and their 
derivatives. The Ui — rLi{r, 9) are hnear in r : Ui — nio(^) + 'n'ii{9)r. In this way, 
by differentiating p.ip two times with respect to r, we obtain a Hnear ordinary 
differential equation for R{r): 

m^R^^^ + (2rn3^ + m2)R^'^^ + (mg^ + 2TO2r + mi)R^^^ 

+ (TO2rr + 2mir)R" + TTlirrR' = (3.2) 

In this expression, by using the form of the Fi{r,9), the coefficients associated 
to the derivatives of different orders of R{r) can be expressed as combinations 
of trigonometric functions and different powers of r. Thus, by separating (j3.2p . 
according to the different hnearly independent trigonometric functions, we ob- 
tain six differential equations for R{r) that have to vanish independently of each 
other: 

Ai(r'*i?(5) _^ 7r3i?(4) _ ^2^(3) _ ;^8^^// _,_ ^8^/) ^ Q (3 3) 

A2{r^R^^^ + 7r3i?(4) _ ^2^0) _ ig^^" _,. ig^') ^ q (3 4) 

Bi{r^R^^'^ + lAr^R^'^'^ + A^r'^R^^'^ + 2ArR" - 2AR') = (3.5) 
B2{r^R^^'^ + lAr^R'^'^'^ + Ad,r'^R^^'^ + 2ArR" - 2AR') = (3.6) 
Cir^ + A3r'')i?("^) + (20Cir'^ + llA^r^^R^'^'' + (l20Cir* + 27 A^^r^^ R'^^^ 

2A0Cir^ + GAgr) i?" + (^I20Cir^ - GAgr) i?' = (3.7) 

(Car^ + A4r'')i?('^) + (20C2r'^ + llA^r^^R^'^'^ + (l20C2r4 + 27 A^r^^R^^^^ 

+ {2A0C2r^ + QA^r^ R" + (l20C2r^ - GA^r^j i?' = (3.8) 

The solutions of these equations are: 
I. ^2+^2 ^0 

R{r) + ttir^ + asr* (3.9) 

r 



( 



2. + 7^ 



3. + C2V 



4. A§ + A| 7^ 



i?(r) = ^ + ^ + a4r2 (3.10) 



^W = ;j + ;3+V (3-11) 



R{r) = — + a4r^ + ae log r (3.12) 
r 
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5. (Ci^a) 

X as , a? , as , / + V A3 + Cir'^ 
R{r) = -7 H - log 



(3.13) 



6. (C2A4) ^ 



„, . 03 ay as . ( \fM + + Cir^ \ 

mr) = = H = log 

r V^4 + C2r2 + C2r2 r ^ 

(3.14) 

We omit \ and the constant terms since they can be absorbed into the angular 
part. 

Theorem 1. A third order integral of motion for a potential of the form 
must have the form and can exist only if one of the following situations 

occurs: 



1. The constants in i2. 9\) satisfy 

Ai=A2 = A3 = A4 = Bi = B2 = Ci = C2 = (3.15) 
Then the radial equations h3.3^) to h3.^) impose no restriction on R{r). 

2. R{r) ^ f fora^O. 

3. R{r) = ar^ for a ^ and Ci = C2 = 0. 

4. i?(r) = 0. 

In cases (ii) and (iv) equations \3.k^) to \3.8\) are satisfied identically for all 
values of the constants in ^L^j) or ^2^ 



Before prooving Theorem 1 let us stress that it gives only necessary con- 
ditions for the existence of the integral F, not sufficient ones. Those will be 
obtained when integrating the equations for the angular part 

Proof, (i) Let us start with Case [i] of the theorem. The condition (3.15) 
implies that (13. 3p to p.Sp are satisfied identically for all R{r). Then (2.5) to 
(2.8) simplify to 

Fi = F3 — 0, F2 — Bq, F4 = + Do 



and we can integrate (2.2)-(2.4) to obtain 



G,{r,9) = -^+m (3.16) 



G2{r,9)=2Bo{R+^) -^S+^+m (3.17) 
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where f3{9) was introduced in (2.14) and ^{9) appears as an integration constant 
(an arbitrary function of 9). Equations (2.2) to (2.4) further imply 

/3 + /3 = 0, ^-3i:>o5' = 0, Bo(S'(^) -45) = (3.18) 

and the integral (2.9) reduces to 

Y - DoLl + BoLM+pD + {gii^: y),Pi} + {92{x, y),P2} (3.19) 

(ii) From now on we assume that at least one of the constants Ai , Bj , Cj 
does not vanish. From p.9p - (|3.14[) we see that the most general form of the 
radial term R{r) is 

/ \ '^l ^2 '^3 9 4 

it(r) = — j- + ^ H ha4r + a^r +a6logr 

a? as ( ^/A + ^/A + Cr"^ 

^ + Cr2 ^ + Cr2 °^ 1^ r 

where (A, C) = (^3,^1) or (^4,(^2). 

Let us show that the "exotic" tems 01,02,6(5,06,07 and og are actually ab- 
sent, i.e. their presence is not allowed by the original determining equations 
(|2.1l) - (|2.4p . We proceed systematically by assuming the contrary. 

(h.l) 01 ^ 0. 

The function R{r) must have the form (3.11) and the only possible nonzero 
constants are Ci and C2 (and as always _Bo and Do). 
From ((221) to (gH), we find 

Fi=0, F2 = Bo, F3^CiCos9 + C2sin9, 
B 1 

Fi = + -(-Ci sine + C2 cos 6*) + Da (3.21) 
r- r 

The term in the compatibility condition (2.10) is 

oi (Ci cos 9 + C2 sin 9)^0 (3.22) 
The condition + C| ^ implies oi = 0. 

(h.2) 01 = 0, 02 7^ 0. 

The function R{r) must have the form p.lOp and the possible nonzero con- 
stants in Y are Bi, i?2, Ci, C2 (in addition to Bq and Dq). 
The coefficient of r° in ((?1U1) is 

O2(BiC0s26'-HB2sin26') = (3.23) 

We can keep 02 7^ only if we impose Bi = B2 = 0. 
The coefficient of r° in (P?T|) is 

02^0^ = (3.24) 
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For J5o = or 5 = 0, the term in in (2.10) implies 

a2(Ci 008 6* + Ca sin 61) = (3.25) 

a contradiction with the assumption + C| ^ (since we aheady have 
Bi = B2 = 0). Thus we have 02 = 0. 

(ii.3) ai — a2 — 0, 05 7^ 0. 

The function R{r) satisfies (j?^ and we impose Al + ^ 0. (PTTi)) and 
(|2.4I) give the exphcit form of the functions Gi and G2 and from the term r^'^ 
in (|2.ip we obtain 

a5(AiCOs36' + A2sin36') = (3.26) 
a contradiction. Hence we have as = 0. 

(n.4) ai = 02 = as = 0, ag 7^ 0. 

The function R{r) must have the form (I3.12p and we request + A^ 0. 
We obtain from (|2.14p and (|2.4p the functions Gi and G2 and (|2.ip in this case 
contains an fog r term with coefficient 

Og (A3 cos 6* + A4 sin 6*) = 

and hence ag = 0. 

(n.5) fli = 02 = as = flg = 0, ag 7^ 0. 

The function R{r) must have the form p.l3p or (I3.14p and we request 
A3C1 7^ or A4C2 ^ 0. Inserting Gi and G2 obtained from ([2TT4)) and ((2^ in 



(|2.1I) we have the term r^^ fog with coefficient 



a| AG"* cos 61 = (3.27) 
Since we impose AC 7^ 0, we find ag ~ 0. 

(ii.6) ai = a2 = as = ag = ag = 0, ar 7^ 0. 
The radial part is 

Rir)^^ + ^=^L= (3.28) 

r VA + Cr^ 

and we must request that AC 7^ 0. From the coefficients of and in (|2.1I) 
we have that: 

arG^ (Bos' - 2a3 A cos 6*) = (3.29) 

ar AG^ (BEOS' - ISaa A cos 6*) = (3.30) 

Then for any values of _Bo, we have 

a3Acos6' ~ 
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Since A we have 03 = 0. In (12.11) we have an term with coefficient 

a7A'*(sin 0^ + 2 cos 6'S') ==0 (3.31) 

and hence S — — 2 cot 9S. Substituting 5 into (12. ip we use the coefficient of r^" 
to obtain 

m = ^^i-h' + &S) (3.32) 
Then the coefficient of gives 

ar^C^ cos 6* = (3.33) 

Since AC 7^ it follows that 07 = 0. 

(iii) So far we have shown that the function R{r) must have the form 

R(r) = — + air'^ (3.34) 
r 

To complete the proof of the theorem we must show that either 03 or 04 must 
vanish. From p.9p - p.l4p we see that the constants in Y that can survive (in 
addition to Dq and Bq) are Ai, A2, A3, A4, Bi, B2. From (|2.ip we obtain the 
coefficient of to be 

04(^1008 36* + A2 sin 30 + yl3 cos 6* + ^4 sin 61) = (3.35) 

and hence for 04 ^ 0, we have Ai — A2 = A^ — A4 — 0. 
The coefficient of yields 

04/3 = (3.36) 

and hence /3 — 0. Taking this into account we return to ()2.10p and find the 
coefficient of r^: 

a3(Bisin26l-B2Cos26l) =0 (3.37) 

For 03 ^ this implies Bi = B2 = and we are back in the generic case where 
the integral is (3.19). 

Thus cither 03 or 04 in (3.34) must vanish and this completes the proof of 
the theorem. □ 

We see that a third order integral Y of (|l.ip with at least one nonzero 
constant Ai, Bj, Cj with i — 1,2, 3, 4, j = 1, 2 can only exist if the radial part of 
the potential is a harmonic oscillator i?(r) = ar^, a Coulomb-Kepler potential 
R{r) = a/r or i? = 0. 



4 Angular term in the potential 

Let us return to the problem of solving the determining equations (|2.ip to (|2.4p 
and concentrate on the angular part, once the radial part is known. We shall 
consider each of the four cases of Theorem 1 separately. 
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4.1 Radial equations satisfied for all R{r) 

According to Theorem 1, the constants in the third order integral Y satisfy 
(3.15) and Y itself is as in (3.19). In this case the functions Gi{r, 9) and G2{r, 0) 
are as in (3.16) and (3.17) where /3(6'), ^{0) and 3(0) satisfy (3.18). 
Two cases must be considered separately: 

1. Bo ^ 

From (3.18) we obtain 

m = iDoS + Co 

^{6) = ^1 cos0 + /32sin6l 

S{9) ^ si 00829 + S2sm2e + So (4.1) 

where ^O: /32, Sq, si, S2 are constants (we can put sq =0). 

Substituing (4.1) in (j2.ip we obtain si = S2 = 0, and we have a purely radial 
potential. In the integral (3.19) we have gi = .92 = and the result is trivial. 
Namely, since L3 is an integral, and L^H are also integrals. In general this 
potential is not superintegrable but first order integrable. 

2) Bo = 0, i^o = 1 

The third equation in (3.18) is satisfied trivially so S{6) in (|4.ip is arbitrary. 
Putting I3{9) and ^(6*) in we obtain 

(/3i cos 61 + (32 fim9yR' ^r{ — S^^'> - 3SS ~ ^qS) 



+ 



We distinguish two subcases: 



((/3i sin 6* - /32 cos 61)5' + 2(/3i cos 9 + P2 sin 61)5') 

(4.2) 



2.a) (3i= 132 = 
We have 

^^25.(3) ^ ;^255 + 4^05" 

This can be integrated to: 

h^S^ = 45^ + 4^05"^ + bS + c (4.3) 

where b and c are integration constants. We can set ^0 = by a suitable change 
of variables S{9) ^ T{9) - ^ and g31) is simplified to 

where t2 and are constants. 

The potential is expressed in terms of the Weierstrass elliptic function: 

V{r^9)=Rir)+'^^^^^ (4.4) 
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for an arbitrary radial part R{r). 

This potential has a third and a second order constant of motion of the form: 

Y ^2Ll + {L3,3h^pi0)} (4.5) 
X ^ Ll + 2h'^p{0) (4.6) 

The second order constant of motion is known as the one-dimensional Lame 
operator. Y and X are algebraicahy related |20) : 

This is a good example of what is called algebraic integrability [121 HO] • Since 
we are looking for hamiltonians with algebraically independent second and third 
order integrals of motion, this result is not a real superintegrable system. 

The classical potential corresponding to (14. 4p is a purely radial potential 
V{r) = R{r). This result is directly obtained from (|4.2I) in the limit ft i— >■ 0. 

2.b) + 
Differentiating (|4.2p two times with respect to r, we obtain: 

and hence 

„/ N a b 
^(^) = - + ^ 

The ^ term can be included in the angular part of the general form of the 
potential. So without loss of generality, we can set b — 0. This case will be 
studied in section 4.2 since the radial part of the potential is of the form of case 
(ii) of Theorem 1 . 

4.2 Potential of the form V(r, 9) = - + with a ^ 

We consider ()2.2p and (|2.4p separate different powers of r and obtain equations 
relating I3{9) in (|2.14|) with the angular part S{0) of the potential: 

/3 + /3 - 2(Ci cos 61 + 6*2 sin d)S + 5(Ci sin 9 ~ C2 cos 9)3 
+ 2(Ci cos 9 + C2 sin 9)S ^ 6a{Bi sin 261 - B2 cos 261) (4.8) 

(Si cos 26* + B2 sin 26* + Ba)S^^^ + 8(-Si sin 261 + B2 cos 26')S' 

+ 4(-5Bi cos 26* - 5B2 sin 26* + Bo)'5' + 16(Bi sin 261 - B2 cos 261)5 

= 3a(-15AiCOs36'- 15^2 sin36' + A3 cos0 + A4 sin0) (4.9) 
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(3^1 sin 36* - 3^2 cos 36* + A3 sin6i - cos 9)8^^^ 

+ (36 Ai cos 36* + 36^2 sin 36* + 4A3 cos 6 + 4^4 sin 9)8 

- (132y4i sin 36* - 132^2 cos 361 - 4A3 sin 9 + AA^ cos 9)8 

- (144Ai cos 361 + 144^2 sin 361 - 16 A3 cos 9 - I6A4 sin 9)8 = Q (4.10) 

Integrating (|2.4p for G2 we obtain the function 

^{9) = ?.Dq8{9) - a(Ci sin 6* - C2 cos 6*) + Co 
From (HH), it follows that S{9) and P{9) verify: 

h^{DQS^^^ + aCi cos6' + aC2 sin6l) - 4^5" + 4a/3 = (4.11) 

^ (— Ci sin 6* + C2 cos 6')S'(3) - 4(Ci cos6' + C2 sin 6*)^ 
+ 6(Ci sin 61 - C2 cos 9)8 + 4(Ci cos 6* + C2 sin 9)8 
+ 6a(-Bi sin 26* + B2 cos 261) ) = 45/3 - 85^ 



- 8(Ci cos 6* + C2 sin 9)8"^ + 12a{Bi cos 26* + B2 sin 26* + Bo)S 

- 12a^(AiCOs36' + A2sin3e' + A3COs6' + A4sine') (4.12) 

{ {Bi COS 261 + B2 sin 29 - Bq)8'^^^ + 8{~Bi sin 26* + B2 cos 261)5 



- 4(5Bi cos 26* + 5B2 sin 261 + Bo)S + 16(Bi sin 261 - B2 cos 26')5 

^ Sah^ {5Ai cos 36* + 5A2 sin 36* + A3 cos 6* + A4 sin 6*) 

+ 12a(3Ai sin 36* - 3^2 cos 36* + ^3 sin 6* - A4 cos 9) 8 

+ 36a(Ai cos 36* + A2 sin 36* + A3 cos 6* + A4 sin 9)8 

+ 2{Bi cos 29 + B2 sin 26* + Bo)55 - 12(Si sin 261 - B2 cos 29)8"^ 

-m{BiCos29 + B2sm29 + Ba)SS (4.13) 
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{ (3Ai sin 36* - 3^2 cos 36* - 3 A3 sin 9 + 3A4 cos 9)8'^^^ 



+ (36Ai cos 36* + 36^2 sin 36* - I2A3 cos 6* - I2A4 sin 9)8 
(-132^1 sin 361 + 132^2 cos36'- 12yl3 sin6' + 12A4 cos 6')5' 

(-144^1 cos 36* - 144^2 sin 36* - 48A3 cos 9 - 48A4 sin 9)8 

- 72 ( Ai cos 30 + A2 sin 36* + A3 cos 6* + yl4 sin 6*) S*^ 

+ (- 120 Ai sin 361 + 120 A2 cos 36* - 4OA3 sin 9 + 40 A4 cos 9) 88 

+ (54Ai cos 36* + 54A2 sin 36* + 6A3 cos 9 + 6 A3 sin 9) S'^ 

+ (6 Ai sin 361 - 6 A2 cos 361 + 2 A3 sin 6*- 2 A4 cos 61)55' (4.14) 

The principal result that fohows from the compatibility of the preceding 
determining equations is that the only potential that satisfies (4.8) to (4.14) is 

+ (4.15) 

r cos"^ 6 

or potentials that can be rotated to (|4.15p . 

The potential (|4.15p is a well known quadratically superintegrable one [31 13] 
with the Coulomb potential as a special case. The third order integral is the 
commutator (or Poisson commutator) of the two second order ones. For more 
recent discussions of the potential (|4.15p see e.g. 13 US- 

4.3 Potential of the form V{r, 6) = ar"^ H — with a ^ 

Here, from Theorem 1, the third order constant of motion is of the form: 

Y = DoLl + BoLsipj + pI) + Bi{L3,pj - pD + B2{L^,PiP2} 

+ {5i,Pi} + {ff2,P2} (4.16) 



As in the preceding case, we obtain from (|2.2p to (|2.4p that the angular part of 
the potential has to be a solution of: 

{El cos 261 + B2 sin 261 + Bq)S'^^'^ + 8{-Bi sin 261 + B2 cos 26')5 

+ A{-5Bi cos 261 - 5B2 sin 26* + ^0)5"+ 16(Bi sin 29 - B2 cos 261)5 = (4.17) 
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for £_{e) = 3DoS + and from (lO) : 

h'^DoS'^^'^ ~12DoSS ~4^oS ^0 (4.18) 

( {Bi cos 261 + B2 sin 29 - Bq)"?'^' + 8{-Bi sin 26 + B2 cos 261)5' 



- 4(5Bi cos 26* + 5B2 sin 261 + Bo)5 + 16(Bi sin 261 - B2 cos 26')S' 

= 2{Bi cos 26* + B2 sin 26* + Bo)SS - 12{Bi sin 26* - B2 cos 26')S'^ 

- 16(BiCOs26' + B2sin26' + Bo)55 (4.19) 

for /3(0) = 0. 

We can solve the compatibility between all those determining equations for 
iS* to obtain the following potential: 



V{r,e) = + 



2 , 2{b + c) + 2{c-b)cos20 



r2 sin^ 26* 



b 

j,2 



ar^ + — + — (4.20) 



or potentials that can be rotated to (|4.20p . 

The potential (|4.20p is also a well known quadratically superintegrable po- 
tential having the harmonic oscillator as a special case (SJ |4| . The third order 
integral can be obtained as a commutator of the second order ones. 

4.4 Potential of the form V{r, 9) = — y- 

From (|2.1[) to (|2.4p . we again obtain the determining equations (4.8) to (4.14) 
with a — Q. The third order constant of motion is in its most general form (j2.9p . 

Solving the determining equations for S{9) and I3{9) (with a = 0), we reob- 
tain special cases of results derived in sections 4.2 and 4.3 (without the Coulomb 
or the harmonic radial parts) . In addition to these known cases we obtain three 
additional ones. The first is 

= (4-21) 
sm 69 

This potential is a special case of the rational three-body Calogero system in 
two dimensions and is already known to be superintegrable '32j. 
The third order constant of motion associated to (I4.2ip is: 

s o 2 r, f —3a;'' + 6x^2/2 + y^'\ { l&xy 
Y =Pi- 3piP2 + 2aipi, — „, , „ — > + a< p2, ■ 



, y2(„32.2+y2)2 j' (_3^2+y2)2 

The potential is obtained both in classical and quantum mechanics. 
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The two others cases occur when the third order integral of motion takes 
the form 

Y = Ci{LIp^} + 2DaLl + + {.92,^2} (4.22) 

i.e. Ai ^ A2 ^ = = Bq = Bi = B2 = C2 ^ 0. In this case, to ([^ 
are reduced to the system for S[9) and (i{0): 

/3 + /3 + Ci (-2 cos 6*^ + 5 sin 6S + 2 cos OS) = (4.23) 
h^DoS^^^ - 4^5* = (4.24) 

h^Cii^- sin 6»S'(3) - 4 cos eS" + 6 sin e^S" + 4 cos 9S^ = 45"^ - SSP - 8C1 cos ^fi-^ 

(4.25) 

with ^{9) = SILOS' + ^0. 

We distinguish two subcases: 

4.4.a) Do = and Ci = 1. 

In this case, we have ^ = or 5 = 0. If we set ^ = and S = T, (|4.23p is 
solved directly for /? 

/3(6l) = ^1 cos 6* + /32 sin 6* - sin 6* + 2 cos Of (4.26) 

Inserting (|4.26p in (|4.25l) we obtain a fourth order ODE for T: 

( sin erf'^) + 4 cos 0r(3' - 6 sin 6*2^ - 4 cos OT)- 12 sin 6irf 
- 4 cos 6irf - 4 (^1 sin 6* - ^2 cos 61) T - 16 cos Of^ 
+ 8 sin6'Tr - 8(^i cos 61 + /32 sin 6*) T = (4.27) 

Under the transformation {d,T{9)) i-> (z,T(z)) where z = tan6' (|427)) be- 
comes: 

h^z{l + z^ fT"" + 4ft2(l + ^2^(1 + 3z2)T"' + [2^2^(13 + 18^2) _ 4^^^ _^ 

- 4r - 12z(l + z^)T')\T" ~ 8(2 + iz^)T'^ - SfiiT' + 4ft^(l + Qz'^)T' = 

(4.28) 

This equation can be integrated twice. The first integral is: 

h^z'^{l + z^fT'" + 2;i2z(l + z2)(l + 3z2)T" - 62^(1 + z2)T'2 + 2{-h^ + h^z"^ 
+ 3^2 z'' - 2/3iz2 + 2/322 ~ 2zT)T' + 27^ - 4^2^ = i^i (4.29) 

for an arbitrary constant of integration Ki . 
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The second integral is: 



^ z(2z^ + i)r - r + /32 



Z2(z2 + 1) 



4z*(z^ + lY{T') 



r\3 



K^z^iz^ + 1)3 

+ z^iz"^ + 1) [4zT + h^iiz"^ + 1) + 4^1^^ - 4/32^] (T')^ - 2z(z2 + 1) [2zT'^ 

- {Af32Z + h^)T - {Kiz - /32h^)]T' - [izT^ + {h^iz^ - 1) + i^iz'^ - 12l32z)T^ 

- 2{Kiz + /32(4/3i + h^)z^ - 4(3lz ~ /32h'^)T - (ft^^az^ - 2l32Kiz + Plh^)] 

(4.30) 

for a second arbitrary constant of integration K2. 
The transformation (z,r(z)) (x, VF(x)): 



2^/xy/l — X 
l-2x ' 



^h^W + {h^ + 4.l3i){l - 2x) 

%\/x\/l — X 



(4.31) 



maps (|4.28p . (|4.29p and (I4.30p to equations contained in a series of papers by 
C. Cosgrove on higher order Painleve equations. Specifically, (|4.28l) is mapped 
into the fourth order equation F-VII [33] (see also [34]). Equation (|4.29l) is 
mapped into the third order differential equation Chazy-I.a of [35] and (|4.30p 
into the second order differential equation of second degree SD-I.a of with 
parameters 



Ci = C4 = C5 = Cg = Cg = 0, C2 = -C3 = 1, C7 = 

ft4 - I6/32 - IQpl - 8K1 



8/3i 



16^12 



cg 



64/14 



ClO 



le/i^Xz - 8(4/?i + ff )Ki + ff{4:/3i + h^f 
256^6 



(4.32) 
(4.33) 
(4.34) 



SD-I.a is the first canonical subcase of the more general equation that Cos- 
grove called the "master Painleve equation", SD-I [36]. SD-I.a is solved by the 
Backlund correspondence 



W{x) 



x-'ix-lf 



APe{Pe - 1){Pg - x) 



P' 

^6 



^^6(^6 - 1) 



-472! 1 



2x 



x{x — 1) 
1 / 2(^-1) 



1 



(l-v/2^)'(l-2P6 



74 

4 



2x{Pg - 1) 



and 



W'{x) = 



x{x — 1) 
4Pe(Pn-l) 



x{x — 1) 



(4.35) 

72(^6 - x) _ 73 (Pg - x) 
2(x - 1)P6 2x(P6 - 1) 
(4.36) 
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where -y/271 can take either sign and 71 , 72 , 73 and 74 are the arbitrary pa- 
rameters that define the sixth Painleve transcendent Pq obtained from the well 
known second order differential equation: 



P" - 



1 



1 



1 



1 

Pe{Pe-l)iPe-x) 



1 



1 



1 



71 



X X — 1 Pq — X 
I2X 73(3^ - 1) lix{x - 1) 

P2 + (Pg _ 1)2 + (pg _ 



P' 

^6 



(4.37) 



The parameters 71,72,73 and 74 are related to the arbitrary constants of 
integration /3i , /32 , Ki and K2 



4c7 = 71 - 72 + 73 - 74 - V 271 + 1, 
4c8 = = (72 + 73)(7i + 74 - V271), 

4c9 = (73 - 72)(7i - 74 - V^Ti + 1) + ^(71 - 72 - 73 + 74 



4cio = ^(73 - 72) (71 + 74 - \/^if + ^(72 + 73)^(71 - 74 



(4.38) 
(4.39) 

(4.40) 

(4.41) 



Only three parameters of (|4.37p are arbitrary in our case. From (|4.39p . we 
see that one of the following relations must hold 



72 = -73, 74 = -71 
From the inverse transformation x z = tan t 



,2 (e\ 



1 



1 



2VI + z2 



(4.42) 



(4.43) 



we obtain two solutions for 5. By taking the derivative of T in (4.31) we obtain 
the quantum potentials 



v(r,e) 



1 U2..rn ^ ±8ti^COs6W{x±)^'ipl+ti^ 

— \nW [x±) , . 2n 

r-' \ 4 sm 9 



(4.44) 



In the limit 1— )■ 0, (|4.29p is reduced to a first order differential equation of 
second degree in T': 

3z^{l + z^)T'^ + 2zTT' -T^ + 2{fiiz^ - (32z)T' + 2P2T + ^ = (4.45) 



(|4.45p is a special case of the more general equation: 

A(z)T'2 + 2B{z)TT' + C{z)T'^ + 2D{z)T' + 2E{z)T + F{z) = (4.46) 
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A number of papers has been devoted to the integration of (|4.46l) . For 
example in ^37j a method is suggested for its integration. 

Special solutions can be obtained under the condition that 



A B D 
B C E 
D E F 



= 



(4.47) 



This condition implies that /?i = and Ki = — 2/3| in ()4.45|) . In this case 
(|4.45l) can be factorized 



Vl+3z^)T' -T + I32) [{z - z^/A + i^)T' - T + j = (4.48) 
We obtain two solutions: 

z3(5 + 3z2 + 2V4 + 3z2)h 



Ti=p2 + a 

h + OL 



(2 + V4 + 3z2)f 
(l + z2)i(2 + V4 + 3z2)! 
z(5 + 3z2 + 2V4 + 322)h 



(4.49) 
(4.50) 



where a is an integration constant. The angular part of the potential is obtained 
by differentiating the preceding results with respect to 9 and the resulting clas- 
sical potentials are 



3q! sec** 



V = 



and 



V = 



7 + 3\/4 + 3 tan2 + cos 26'(l + x/I+Tta? 



r-2 tani 6IV4 + 3 tan2 {2 + V4 + 3 tan2 9^ ' + 3 tan2 + 2\/4 + 3 tan2 9 



a sec 3 



(4.51) 



47 + 17\/4 + 3tan2 6l + 18 cot2 9{2 + V4 + 3tan2 9) + 3 tan2 61(5 + x/ITsIm? 



2r2V4 + 3tan2 6l(^2 + \/4 + 3 tan2 9^ ' ^5 + 3 tan2 9 + 2V4 + 3 tan2 9 



(4 



If condition (14.471) is not satisfied, the general solution of (|4.46p is related to 
the general solution of the equation 



dw M{z)w'^ + N{z)w'^ + P{z)i 



dz 



where 



T{z) 



+Q{z) 



w + mw + n 
pw 



(4.53) 



(4.54) 
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and Af , N, P and Q are complicated algebraic expressions depending on z and 

A B D 



2(BD ~ AE) 

m = =, n 

ylVS2 - AC 



-1 



A{B^ - AC) 



B C E 
DEE 



P 



2VS2 - AC 
A 



4.4b) Do = liDo^ 0). 

Equation (4.24) can be integrated to (4.3) so the angular part of the potential 
is expressed in terms of the Weierstrass elliptic function S{9) = h?p{9). The 
corresponding potential is 



m 



vir,e) 

I - 

From the compatibility of (4.23) and (4.25) we obtain 
-Ci 



(4.55) 



(6^ + 5c(n2 _ 2^0)) cos e - 2(36 - ft-* 
^ + 6(;i^+eo)sin6lp+(-2(9c4 



(4.56) 



4(/i2 + Co)p^+46p + 6c 
+ 2^^) cos 6'p2 - 8 (ft^ + ^o) cos ( 

+ Co))cos6l + (36 + 2ft2(fi2 +^o))sin6'p)p 
The potential ()4.55p thus allows two third order integrals of motion (see 

Km ) 

Fi - 2^3 + 1^3, 3nV(^)}, (4.57) 
Y2 = {L^, pi} + cos 6* + (2 cos 6'p-/3) sin 6*, pi} 

+ {P sin 6* - (2 cos Op - /3) cos 9, P2} (4.58) 

where /3 is the expression in (j4.56p . 

The integral (|4.57l) coincides with (|4.5p for the more general potential (|4.4p . 
As noted above (|4.4p is not really superintegrable because of relation (|4.7I) . The 
potential (j4.55p is superintegrable since I2 (|4.58p . the second order integral of 
motion X (|4.6I) and the Hamiltonian H are algebraically independent. 

In the classical limit hi—>-0, the system reduces to free motion. 



5 Summary and Conclusion 

The main results of this study can be summed up in two theorems. 

Theorem 2. In classical mechanics in the Euclidean plane precisely 4 classes of 
Hamiltonian systems separating in polar coordinates and allowing a third order 
integral of motion exist. The corresponding potentials are (4-15), (4-20), \4-21^ 
and 

Vir^e)^^ (5.1) 

where T{z) satisfies equation \4-45^ for z — tan0. The third order integral of 
motion is 

Y = {L^,pi}+{;9cos6i+(2cos6if-/3) sm9,pi}+{/3 sin^— (2 cos OT—p) cos 9,p2] 

(5.2) 
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The potential (4.15) and (4.20) are quadratically superintegrable and well- 
known. The third order integral is functionally dependent on the quadratic ones. 
The potential (I4.2ip is the three body Calogero system with no central term. 
Thus the genuinely new superintegrable classical potential is (5.1). We have not 
obtained the general solution of equation (|4.45p , but particular solutions led to 
the superintegrable potentials (4.51) and (4.52). 

Theorem 3. In quantum mechanics the superintegrable systems correspond to 
the three known potentials (4.15), (4-20) and \4-21^ plus two new ones. One 
new one is given by (4.44) where W and W' are expressed in terms of the sixth 
Painleve transcendent Pg in (4-35) and (4.36). The other new one is given by 
the Weierstrass elliptic functions p{9) in (4.55). 

The Painleve transcendents were first introduced in a study of movable sin- 
gularities of second order nonlinear ordinary differential equations. They play an 
important role in the theory of classical infinite dimensional integrable systems. 

The transcendent Pq that was obtained as a superintegrable quantum po- 
tential in this article depends on 3 free parameters (see (I4.42p ). The Painleve 
trancendents Pi, P2 and P4 have already appeared for potentials separable in 
cartesian coordinates |22j . A remarkable relation between quantum superinte- 
grability and supersymmetry in quantum mechanics was discovered and used 
to solve the Schrodinger equation with potentials expressed in term of Painleve 
transcendents p5]-[30). 
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